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A method is outlined for calculating the reflection coefficient from a horizontally 
stratified ionized medium. The profiles of electron density and the collision frequencies 
are both taken to be exponential functions. The d-c magnetic field is taken to be horizontal 
and transverse to the direction of propagation. The specific results described are applicable 
to the oblique reflection of VLF radio waves in the D layer of the ionosphere for propagation 
along the magnetic equator. It is confirmed that the reflection coefficient is nonreciprocal 
in both amplitude and phase. For a wide range of the parameters, the magnitude of the 
reflection coefficient is greater for west- to-east propagation than for east- to-west propagation. 



1. Introduction 

The lower ionosphere is primarily responsible for 
the propagation of VLF radio waves to great dis- 
tances. In theoretical treatments of this problem 
it is often assumed that the lower edge of the iono- 
sphere may be represented by a sharply bounded 
and homogeneous ionized medium. Actually, such 
a model was used by G. N. Watson [1919] over 40 
years ago. Applications and refinements of such 
a model have been discussed frequently in the 
recent literature [e.g., Budden, 1962; Wait, 1962; 
Johler, 1963]. One such refinement is the inclusion 
of the earth's magnetic field in the analysis. This, 
of course, renders the medium anisotropic. If the 
vertical inhomogeneity (or horizontal stratification) 
of the ionosphere is also considered simultaneously, 
the situation becomes very complicated indeed. An 
extensive study of analytical methods to treat such 
problems has been carried out by K. G. Budden and 
his colleagues at Cambridge University. Much of 
this work is summarized in a monumental text 
[Budden, 1961] which will be the standard work on 
the subject for some time. Budden makes extensive 
use of "full wave" methods which may be described 
as a frontal assault on the differential equations 
satisfied by the field components in the medium. 

In this paper we shall consider a special case of 
a horizontally stratified and anisotropic ionosphere. 
Specifically, the earth's magnetic field is assumed 
to be purely transverse to the direction of propagation. 
Strictly speaking, this is applicable only to the 
situation when the path of propagation is along the 
magnetic equator. However, the characteristics 
in this special case prevail at other latitudes if the 
transverse component of the field is appreciable. 
At least this is borne out by a numerical study of 
the sharply bounded ionosphere for an arbitrary 
magnetic dip angle [Johler, 1961]. In any case, 
the resulting simplicity of the differential equations 



for the limiting case of a purely transverse magnetic 
field encourages one to consider this situation in 
more detail. In particular, it is desirable to investi- 
gate the influence of gradient of both the electron 
density and collision frequency. In much of the 
previous work on this subject the collision frequency 
has been assumed constant. 

In a study of the recent literature [e.g., Belrose, 
1963] it is found that both the electron density N(z) 
and the collision frequency v(z) vary approximately 
in an exponential manner with height z. For ex- 
ample, in the undisturbed daytime ionosphere we 
may assume that 



and 



N(z)=N exp(bz), 
v{z)=v Q exp (—az), 



(1) 
(2) 



where a and b are positive constants and z is some 
specified level in the ionosphere. From a study of 
the experimental data [Belrose, 1963], it appears that, 
if the reference level is 70 km above the earth's sur- 



face, iVo^-lO 2 electrons/cm 3 and 



*v 



The 



gradient parameters are then expected to be given 
approximately by 6^0.15 km -1 (±0.1) and a ^0.15 
km -1 (±0.02). The quoted values of these con- 
stants must be considered tentative and certainly 
subject to change. Furthermore, it must be under- 
stood that significant departures from the exponential 
shape are to be expected under disturbed conditions. 

2. Formulation 

The situation is shown explicitly in figure la. A 
vertically polarized plane wave is incident at angle 6 
on to a horizontally stratified ionosphere. The z 
axis is taken to be positive in the upward direction. 
At the reference level z=0, the electron density and 
the collision frequency have values designated by No 
and v , respectively. The "scale height' ' which is 
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Figure la. Illustrating oblique reflection of VLF radio 
waves from an inhomogeneous (horizontally stratified) 
ionosphere with a transverse d-c magnetic field. 

equal to 1/b or 1/a, as indicated in figure lb is of the 
order of 6 km for both of these profiles. As men- 
tioned above, these are typical of the daytime D 
layer for both the N and v profiles. 

The lower ionosphere, which is idealized here as a 
stratified ionized medium, may be regarded as an 
electron plasma. The (angular) electron plasma 
frequency co is thus given by 



a^=3.18X10 9 XiV, 



(3) 



where N is the electron density in electrons per 
cubic centimeter and co has dimensions of radians 
per second. The dielectric properties of such a 
(cold) plasma may be described in terms of a tensor 
dielectric constant (e). Thus, the displacement 

vector D and the electric field E are connected by 



D=(e)E. 



(4) 



Choosing the d-c magnetic field to be along the axial 
direction, 1 the tensor has the form 



(e) = 



o 



r 
i 





" 




(5) 



where harmonic time dependence according to exp 
(+icot) is assumed. The coefficients e, e" ', and q 
are given by [Wait, 1962] 
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(8) 



1 In what follows, the x axis is taken to be in the axial direction. Thus eq (4) 
is given by 

/Dy\ /Ey 
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N(z) = N exp(bz) 
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Figure lb. Sketch of the profiles of electron density N(z) 
and the collision frequency v(z). 



where e =8.85X10 -12 is the dielectric constant of 
free space and oo T is the (angular) gyro frequency. 

Maxwell's equations, which are applicable to this 
situation, are 



i(e)ooE=Q,m:\ H, 



iij,coH= curl E, 



(9) 



(10) 



where fx is the magnetic permeability of the plasma. 
The first equation above may be written 



where 
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(14) 



For convenience in what follows, we now choose 
the gyro axis to be in the x direction (i.e., the d-c 
magnetic field is taken to be parallel to the x direc- 
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tion). As a result d/dx=0 since the incident wave 
is specified to be transverse to the d-c magnetic 
field. It is now a straightforward matter to show 
that 

d2 ' * ' F ■ H x =0, 



W ^ 2 mm) 

( d 2 d 2 fcV'X 



E x =0, 



(15) 



(16) 



where &=(e 0) u )^a;=27r/A and M=/z /ju. Here X is 
the free space wavelength and fi is the magnetic 
permeability of free space. 

From Maxwell's equations, the other field com- 
ponents may be obtained from H x and E x . For 
example, 

te Q u)E v =M- ±iK 



U GoE 2 =iK 



—i/j, a)H y 



bz ' 
dH x 



-M 



dy 
dy 



dz 



i/d {) ojH 2 =—M-^' 



(17) 
(18) 

(19) 
(20) 



It is immediately evident that the general prob- 
lem splits into two parts. The total field may be 
regarded as the superposition of two partial fields; 

one is characterized by E=(E X , 0, 0) and the other 

by H=(H X , 0, 0). These may be called TE (trans- 
verse electric) and TM (transverse magnetic) waves, 
respectively. For a vertically polarized incident 
wave, the first partial field (i.e., the TE waves) is 
not excited, so further attention is restricted to the 
TM waves. 2 

3. Preliminary Problem 

As a simple preliminary problem we shall consider 
the ionosphere to be sharply bounded and homo- 
geneous. Thus, for £>0, the dielectric properties 
are to be characterized by the tensor (e) which does 
not vary with z. The region 2<0 is free space. 

The incident wave is of the form 



H' x nG =h e- ikCz e- ikSy , 



(21) 



where h is a constant and where C=(l — S 2 ) 1A = cos 6 
in terms of the angle of incidence d. Then the 
reflected wave must be of the form 

H?«==hoe +ikCz e- ikSy R, (22) 

where R, by definition, is the reflection coefficient. 



2 Actually, the TE waves in this case are not influenced by the d-c magnetic 
field and the medium is effectively isotropic. 



For the region 2>0, 

U x =j(z) exp (-ikSy), 
where, by virtue of (15), j(z) satisfies 



L° z \MM /J 



f(2) = 0. 



Therefore, 

H x =h Q Texp[-ikl-^- 
L LMM 



s 



3 J exp (— % 



(23) 
(24) 

■Sy), 

(25) 



where T, by definition, is a transmission coefficient. 

The unknown functions R and T arc 4 found by 
applying the boundary conditions which require the 
continuity of E v and II X at z = (). This readily leads 
to 

T= 2C 



and 



where 
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where 
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1, which is the usual case 
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The latter result was given by Barber and Crombie 
[1959] who employed a somewhat more involved 
derivation. 

4. The General Problem 

We shall now return to our originally stated 
problem. The continuous profiles of N(z) and v(z) 
are replaced by a very large, but finite, number 
of steps. In other words, the in homogeneous medium 
is replaced by a stack of thin homogeneous layers. 
For purposes of discussion, there shall be P such 
layers while a typical layer is the pth. layer. Tims, p 
ranges from 1 to P through integral values. Some- 
where at a sufficiently negative value of z, the 
medium may be regarded as free space. This level 
is denoted z= — z . 

Within the pth layer, the fields may be regarded 
as a superposition of upgoing and downgoing waves. 
These are characterized by the functions exp 
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(—ifipkz) and exp (+ifi p kz), respectively, where 

P,= -^S 2 T- (30) 

LM P M P J 

The characteristic impedances associated with these 
wave types are 

£?=-S> (3i) 



and 



Ey 

Kj,=jp:> 



(32) 



where the + sign denotes upgoing and the — sign 
denotes downgoing waves. From (17) and (18), it is 
seen that 

K+ = Vo (M p p p +iK p S), (33) 

and 

K;=vo(M p p p -iK p S), (34) 

where 

Vo= (fJLol e ) y2 = 120Tr ohms. 

The problem may now be solved by an application 
of nonuniform transmission line theory [Schelkunoff , 
1943; Wait, 1962]. Thus, the reflection coefficient 
which is referred to the lower edge of the bottom 
slab is given by 

K£-Z x C-A 



Bo 



''Ko+Z x tf+A' 



(35) 



where A=Z 1 /r] and where Z x is the input impedance 
at the bottom of layer number 1. Now, Z x may be 
expressed in terms of Z 2 which, in turn, may be 
expressed in terms of Z 3 . The process is continued 
until the topmost layer is reached where Z P =K P is 
assumed known. The details of this derivation are 
given elsewhere [Wait, 1962]. 

The required number of layers is best determined 
by studying the stability of the solution as the num- 
ber is increased. Because of the relatively long 
wavelength involved and because of the finite losses 
in the medium, the solution converges nicely as the 
number of layers is increased. For the cases dis- 
cussed here, the step size was never greater than 0.1 
km and this appeared to be smaller than necessary. 

5. Presentation of Results 

The final results of the numerical calculations are 
presented in such a fashion that the phase of the 
reflection coefficient R is referred to the level z=0. 
Thus, by definition, 



R=[R ex V (i2Wz )] Zo - 



(36) 



Physically, this means that the observer at £= — z Q is 
sufficiently far below the ionosphere that the medium 
may be regarded as free space. In practice z Q is 
chosen to be large enough that the phase of R does 
not vary with further changes in z . For the case 
here, z was of the order of 40 km. This particular 



normalization of the phase has been used on previous 
occasions [Wait and Walters, 1963]. 

Following the usage in previous papers [Wait and 
Walters, 1963], the quantity o) r =o)l(0)/v(0) is speci- 
fied. In particular, co/co r =K at 15 kc/s, or oo r = 
6-7rX 10 4 sec" 1 . The "effective" conductivity <r e of the 
medium at this level 2=0, is then given by a e = 
e co r ~ 1 .7 X 10" 6 mhos/meter. With exponential-type 
profiles the fixing of the parameter of oo r is not an 
essential restriction. It is a simple matter to shift 
the reference level from z=0 to any other value if 
desired. 

The parameters of the problem are thus X, C, b, a, 
and cc T /v . In order to display the relative influence 
of these quantities, it is desirable to plot the ampli- 
tude and phase of R as a function of co T /v from —3 
to +3 for a range of values of X, C, b, and a. It 
should be noted that X is in km, C is dimensionless, 
while b and a have dimensions of km -1 . Conse- 
quently, the scale length in the present problem is 
the kilometer. By changing this scale, the results 
may also have significance at higher frequencies. 

In figures 2a and 2b the magnitude of the reflection 
coefficient \R\ and the phase of R are plotted as a 
function of a) T /v Q . Negative values of the abscissa 
correspond to propagation from west to east along 
the magnetic equator. The cosine angle of incidence 
is fixed at 0.1. Thus, the angle is highly oblique, 
being only 5.7° from grazing. For long-distance propa- 
gation of VLI radio waves, such highly oblique con- 
ditions prevail. 3 For the curves in figures 2a and 2b, 
the collision profile is chosen so that the collision 
parameter a=0.15 km -1 and the wavelength X=15 
km correspond to a frequency of 20 kc/s. For these 
curves, the electron density parameter b takes the 
values 0.1, 0.2, and 0.5 km -1 . It is evident that for 
oo T /v =0, the steep gradient of electron density is 
associated with maximum amplitudes of reflection. 
However, when co T /v is finite, this may no longer be 
the case. In fact, the asymmetry of the curves about 
lo t /vo=0 is a measure of nonreciprocity in the reflec- 
tion process. As indicated, the reflection coefficient 
for propagation from west to east is greater than for 
propagation from east to west. There is also some 
nonreciprocity in the phase curves but it is not great. 

In figures 3a and 3b, a set of curves show the in- 
fluence of varying the collision frequency parameter 
while keeping the electron density parameter fixed 
at 6 = 0.15 km -1 . For these curves, as before, (7=0.1 
and X = 15 km. It is evident that the steeper gradient 
of the collision frequency corresponds to larger 
reflection coefficients. 

In figures 4a and 4b, \R\ and the phase of R are 
shown as a function of o) T /v for various frequencies 
in the range from 6 to 60 kc/s. For these curves, 
(7=0.1, 6 = 0.15, and a=0.15. The tendency is for 
the reflection coefficient to be diminished at the 
higher frequencies. In this case, the medium is 
acting like a good absorber rather than a reflector. 



3 In fact, the attenuation of the dominant mode in the earth-ionosphere wave- 
guide at VLF is approximately proportional to 1 — | R\ for highly oblique incidence 
[Wait, 1962J. 
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Figures 2a and 2b. Reflection coefficient curves illustrating 
(he dependence on electron density profile. 

In figures 5a and 5b, \R\ and the phase of R are 
plotted for different values of the angle of incidence. 
For these curves, X=15 km, 6 = 0.15 km -1 , and 
a=0.15 km -1 . In general, it may be seen that the 
reflection coefficient is diminished for the steeper 
angle of incidence. It is rather interesting to note 
that the asymmetry (or nonreciprocity) in the phase 
curves is more pronounced at the steeper angles. 

Finally, in figures 6a and 6b, \R\ and the phase of 
R are plotted for various values of the collision fre- 
quency parameter a for \= 15 km and (7=0.1 . These 
curves differ from figures 3a and 3b in that here the 
parameter fi=b J ra is fixed rather than just a. In 
other words, the profile of N/v as a function of z is 
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Figures 3a and 3b. Reflection coefficient curves illustrating 
the dependence on collision frequency profile. 

fixed while the gradient of v is changed. In the iso- 
tropic case, where co r =0, it is interesting to note 
that R is determined only by the gradient of N/v. 
However, for a finite gyrofrequency, the situation 
is changed significantly. In general, the nonreci- 
procity is accentuated when a is diminished. For 
example, if v were assumed to be a constant, the 
dependence of the gyrofrequency is much greater 
than for a collision frequency which varies with 
height. In much of the earlier work [e.g., Budden, 
1 955] on full wave solutions in ionospheric radio waves, 
it is often assumed that v can be regarded as a con- 
stant. Clearly, such an assumption may lead to very 
misleading results. 
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Figures 4a and 4b. Reflection coefficient carves illustrating 
the dependence on wavelength. 



6. Discussion and Concluding Remarks 

The numerical results given here should provide 
some insight into the nature of reflection from an 
inhomogeneous ionized medium. The nature of the 
dependence on N, v, and u T is quite complicated. 
Nevertheless, it appears that the sharper gradients 
of electron density are usually associated with higher 
reflection coefficients. The dependence on the colli- 
sion frequency profile is not so clear-cut. 

In nearly every case it may be seen that the pres- 
ence of the transverse magnetic field is to cause the 
reflection coefficient |jR| to be nonreciprocal. Fur- 
thermore, for a wide range of the parameters, \R\ is 
greater for west-to-east propagation than for east- 
to-west propagation. This is in accord with experi- 
mental data of Round et al. [1925] who observed 
that, for propagation over distances of the order of 
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Figures 5a and 5b. Reflection coefficient curves illustrating 
the dependence on angle of incidence. 

6000 km, signals from VLF transmitters to the west 
are received more strongly than from those to the east. 
This observation has also been confirmed by Crom- 
bie [1958] in a series of field strength measurements 
in New Zealand and by Taylor [1960], who analyzed 
the waveforms of atmospherics. For some inex- 
plicable reason, Budden [1955] deduces from a full 
wave solution that the directional dependence is just 
opposite to this. Although his model is not the same 
as the one considered here, it is difficult for this 
writer to accept the validity of his results in this 
regard. However, it is possible that, because of the 
complexity of the various phenomena, a reversed 
trend may emerge for certain special conditions, par- 



100 



OJ 



IRI 0.6 



0.4 



0.3 







1 1 

X = 


1 

15km 







c = 


0.1 _ 




^5 N. 


$-- 


b + a = 0.3 km" 




-Si^ >A 








0.2 "^ ^< 




-^?o 


- 










N& 








\ \ \ 




1 1 


1 1 


i\ N 



-130° 



-140° 



in 
< 

x 

Q_ 



-150° 



-160° 



-3 -2 




"tAo 



1 




1 1 


,««• 


i J-- rsi 


- o = o.zJSL- 




-^5^^^"^ 




^"^^£25 




*y ' 


X = 15km 




\p 


yS o/ 




C = O.l 
$ = 0.3 km" 1 

1 1 




1 1 ^ 



oj t /i/ 



Figures 6a and 6b. Reflection coefficient curves illustrating 
the dependence on collision frequency profile when the N/*> 
-profile is fixed. 



ticularly for the nighttime ionosphere [Rhoads et al., 
1963]. It is also worth mentioning that Budden 
[1955] has some qualms concerning the accuracy of 
his numerical data at small values of C. 
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